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Abstract. Antiferromagnetic Heisenberg spin chains with various spin values {S — 1/2, 1,3/2,2,5/2) are 
studied numerically with the quantum Monte Carlo method. Effective spin S chains are realized by ferro- 
magnetically coupling n = 2S antiferromagnetic spin chains with 5* = 1/2. The temperature dependence 
of the uniform susceptibility, the staggered susceptibility, and the static structure factor peak intensity 
are computed down to very low temperatures, T/J ~ 0.01. The correlation length at each temperature is 
deduced from numerical measurements of the instantaneous spin-spin correlation function. At high tem- 
peratures, very good agreement with exact results for the classical spin chain is obtained independent of the 
value of S. For the S=2 chain which has a gap A, the correlation length and the uniform susceptibility in 
the temperature range A < T < J are well predicted by the semi-classical theory of Damle and Sachdev. 
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1 Introduction 

For many years, low-dimensional quantum magnets have 
drawn much interest from both the theoretical and ex- 
, perimental condensed matter physics communities. Pow- 
' erful experimental tools like neutron scattering have eluci- 
] dated the basic physics of two-dimensional (2D) antiferro- 
. magnets, such as K2NiF4 (S'=l)[l], and one-dimensional 
; (ID) systems, such as (CD3)4NMnCl3(TMMC) (5=5/2) 
[2,3]. The discovery of high-temperature superconductiv- 
ity in 1986 sparked renewed interest in this field, since the 
[ parent compounds of these superconductors provide very 
I good realizations of 2D spin-1/2 square-lattice quantum 
Heisenberg antiferromagnets[4,5]. Recently, quantum spin 
, ladders and spin chains have also attracted much atten- 
■ tion, since new copper oxides with such structures have 
' become available; these systems are of intrinsic interest 
and they also allow one to compare experimental results 
with the relevant quantum field theories and numerical 
simulations [6]. 

We report here a Monte Carlo study of antiferromag- 
netic spin chains as a function of spin quantum number 
S and temperature T. Recent advances in computer tech- 
nology as well as the development of a powerful loop clus- 
ter algorithm [7] have made it possible to perform a very 
detailed study of various experimentally relevant ther- 
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modynamic quantities. In this paper, we compute the S 
and T dependence of the uniform susceptibility, Xu(S', T); 
the spin-spin correlation length, ^{S,T); the staggered 
susceptibility, Xs(<S', T); and the static structure factor at 
q — TT, 6*^(5, T). There exists a large number of quasi-lD 
magnetic systems, such as Sr2Cu03 (5=1/2) [8], copper 
benzoate (5=1/2) [9], Ni(C2H8N2)2N02C104 (5=1) [10], 
Y2BaNi05 (5=1) [11], CsVCla (5=3/2) [12], (CioH8N2)MnCl3 
(5=2) [13], and TMMC (5=5/2) [2]. These magnets all 
exhibit nearly ideal ID behavior over a considerable range 
of temperature. Therefore, our study facilitates compar- 
isons among experimental results, numerical simulations, 
and theories for such materials. 

The Hamiltonian for the nearest -neighbor Heisenberg 
spin chain is 

H = J^S, -S^+i, (1) 

i 

where J is positive for an antifcrromagnet. We use units in 
which h = ks = gjJiB = 1- The study of the Hamiltonian 
Eq. (1) has a long history that dates back to the remark- 
able exact solution found by Bethe in 1931 for 5=1/2[14]. 
He found the ground-state eigenfunction for this system 
and showed that there is no long-range order at T = 0. 
In 1962, des Cloizeaux and Pearson [15] derived the ex- 
act dispersion relation of the lowest-lying excited states 
at T = 0. Luther and Peschel[16] showed that the spec- 
trum is gapless without magnetic ordering in the ground 
state, and that the spin correlation function decays al- 
gebraically with distance. Haldane[17] eventually conjee- 
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tured that half -odd-intcgcr spin chains would behave qual- 
itatively like a 5=1/2 chain, while for integer spin chains 
the zero-temperature spin correlations would decay expo- 
nentially with distance due to the presence of an energy 
gap, A, in the excitation spectrum. Both numerical and 
experimental confirmation of this conjecture followed [18, 19] . 
The low-energy and low-temperature properties of inte- 
ger spin chains are well described by the ID quantum 
0(3) nonlinear a model without any topological term [17]. 
When conformal field theory techniques were applied to 
ID quantum systems, it was shown that the integrable 
S=l/2 Hamiltonian, at low energies, is equivalent to the 
SU(2) Wcss-Zumino-Novikov-Witten (WZNW) model with 
a topological coupling constant k — 1 [20]. Moreover, 
all half-odd-integer spin-S Heisenberg models were pre- 
dicted to be equivalent to the k = 1 WZNW model, inde- 
pendent of S'[21,22]. 

In the limit of classical spins, 5 — > oo, the Hamiltonian 
Eq. (1) was solved exactly by Fisher in 1964 [23]. In order 
to accommodate the limit 6*^00, Eq. (1) is conveniently 
rewritten in terms of the unit vector Sj = Sj/ y/S{S + 1), 
thus introducing the energy scale JS{S + 1) in place of 
J. The results for the correlation length and the uniform 
susceptibility per spin are, respectively. 
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and 



SiS+l)l + uiS,T) 

Xu(^> 00,T) = — TTTTFT' (3) 



3T 1 - u{S, T) ' 
where a is the lattice constant and u{S, T) is given by 

'JS{S+l) 



u{S, T) = coth 



T 

JS{S+1)' 



For the experimental system TMMC {S = 5/2), the cor- 
relation lengths in the Heisenberg regime obtained from 
neutron scattering could be explained very well by Eq. 
(2), without any adjustable parameters, since J had been 
determined independently [2]. 

In Sec. 2, we give a brief description of our Monte 
Carlo method. Uniform susceptibility data are shown and 
discussed in Sec. 3, while staggered quantities such as the 
correlation lengths, staggered susceptibilities, and static 
structure factor peak intensities are presented in Sec. 4. 
The properties of integer spin chains in the context of 
recent field theory results are discussed in Sec. 5. A com- 
parison between the results for the spin-5 chain and the 
n-chain 5=1/2 ladder is also made in Sec. 5. 



2 Quantum Monte Carlo 

We have carried out quantum Monte Carlo simulations 

on large lattices utilizing the loop cluster algorithm[7]. In 
order to realize chains with spin S > 1/2, an n-chain 
spin-1/2 ladder with an infinitely strong ferromagnetic 
inter-chain coupling is mapped to a 5 = n/2 chain[24]. 
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Fig. 1. The uniform susceptibility per spin is shown as 
a function of T/[JS{S + 1)]. The dashed line is a plot of 
Fisher's result for the classical spin system, Eq. (3). The 
dot dashed lines are fits for 5 = 1 and 5 = 2 to Eq. (4). 
The solid lines for 5 = 1/2 and 5 = 3/2 are the WZNW 
nonlinear a model expression, Eq. (5). 



The same algorithm employed to study spin ladders is 
used with minor modifications[25]. The lengths and Trot- 
ter numbers of the chains are chosen so as to minimize any 
finite-size and lattice-spacing effects. The chain length is 
kept at least 10 times larger than the calculated correla- 
tion length. Spin states are updated about 10^ times to 
reach equilibrium and then measured 10^ times. We are 
able to simulate as many as 100 million lattice points on 
our workstation. 



3 Uniform Susceptibility 

The uniform susceptibility is shown in Fig. 1 as a function 
oiT/[JS{S+l)]. The dashed line is a plot of Fisher's result 
for the classical spin chain, Eq. (3). At high temperatures, 
T/[JS{S+1)] > 1, the agreement among the results for all 
the spin values and Eq. (3) is very good. As predicted by 
Haldane[17], at low temperatures the uniform susceptibili- 
ties of integer spin chains behave markedly differently from 
those of half-odd-integer spin chains. An exponential ac- 
tivation due to an energy gap is observed for integer S. By 
fitting the uniform susceptibility to the low-temperature 
expression [26,27] 
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where A is the Haldane gap and v is the spin-wave veloc- 
ity, we extract the values As=\/J = 0.40(1), vs^i/Ja = 
2.5(1), As=2/J = 0.090(5), and vs=2/Ja = 4.50(33). The 
uniform susceptibility data are fit only for T < A/2, and 
the results of this analysis agree with the values deduced 
in previous studies[18,28,29,30]. 

For the 5 = 1/2 chain, the theoretical low-temperature 
result for the WZNW nonlinear a model [31], 



Xu(l/2,T) 
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is plotted as a solid line. Here the spin-wave velocity 
vs=i/2/Ja = 7r/2, and To/ J « 1.8. The hypothesis that all 
half-odd integer Heisenberg chains arc in the same uni- 
versality class as the 5=1/2 chain, namely the A: = 1 
WZNW model, is now generally accepted. However, this 
claim was initially considered to be controversial, since in- 
tegrable spin-S Haniiltonians are equivalent to the WZNW 
model with k = 25" [32]. A recent numerical study by Hall- 
berg et al. [33] gave strong evidence that fc = 1 for the 
WZNW model of the 5 = 3/2 chain, thus supporting the 
above hypothesis. We give here additional evidence in fa- 
vor of this claim. The asymptotic value of the uniform sus- 
ceptibility for the k = 2S WZNW model is different from 
that of Eq. (5). Specifically, at T ^ 0, Xu{S,T) -> S/wvs 
rather than Xu('S', 0) = l/27rus. If we use the calculated 
value U5=3/2/Ja = 3.87 for the spin-wave velocity[33], 
we obtain (3/2,0) ~ 0.12 for the k = 2S model and 
JXu (3/2,0) » 0.04 for the = 1 model. Therefore, as is 
evident from Fig. 1, our data for 5 = 3/2 show clearly 
that the 5 = 3/2 Heisenberg chain is equivalent to the 
WZNW model with fc = 1. We have, in fact, fitted our 
5 = 3/2 chain data for Xu with Eq. (5). The spin-wave 
velocity is held fixed at the value given above and only Tq 
is adjusted. As may be seen in Fig. 1 the fit is very good at 
low temperatures, as expected. For the 5 = 5/2 chain, we 
are not able to obtain numerical data down to low enough 
temperatures to yield a meaningful test of Eq. (5). 



4 Correlation Length, Staggered 
Susceptibility, and Static Structure Factor 
Peak Intensity 

In order to deduce the spin-spin correlation length, ^(5, T), 
the instantaneous spin-spin correlation function, C(r), is 
computed and fitted to the asymptotic form 



C{r) 



exp(-r/0 



(6) 



Equation (6) is equivalent to ID and 2D Ornstein-Zernike 
(OZ) forms when A = and A = 1/2, respectively. Only 
data with r > 3^ are included in the fit to ensure that the 
asymptotic behavior is probed. For half odd integer spin 
chains, the ID OZ form is found to work very well over the 
entire temperature range. The spin correlations for integer 




T/[JS(S+1)] 

Fig. 2. The spin-spin correlation length deduced from fit- 
ting the computed correlation hmction to the asymptotic 
form, Eq. (6) is plotted as a function of temperature. Also 
shown are the correlation lengths of TMMC from Ref . [2] . 
The classical spin system result, Eq. (2), is shown as a 
dashed line. The sohd line is Eq. (7) for the 5=1/2 
chain. 



spin chains, however, are found to go through a crossover 
with decreasing temperature, from the ID to the 2D OZ 
form at T « A/2. Specifically, the 2D OZ form gives a 
better description of the computed spin correlations at low 
temperatures. This behavior also was observed previously 
for spin ladders of even width[25]. The low-temperature 
correlation function for the 5 = 1 chain is known to be 
proportional to the modified Bessel function Ko{r / £_)[29]; 
Kq is also used to fit our integer spin chain data, and the 
results are consistent with those extracted using Eq. (6). 

In Fig. 2, the numerical correlation length is shown to- 
gether with £,{T) for TMMC obtained from neutron scat- 
tering experiments [2]. The result for the classical spin sys- 
tem, Eq. (2), is shown as the dashed line. The agreement 
among the spin correlation lengths for the different quan- 
tum spin chains, and with the classical curve, is very 
good for temperatures as low as T/[JS{S + 1)] « 0.2, 
which corresponds to ^/a « 4. As the temperature is low- 
ered further, the data begin to deviate from the classi- 
cal curve. Note, that for 5 > 1, as the spin value is in- 
creased, the agreement between the quantum and classi- 
cal results persists down to progressively lower tempera- 
tures. Because of the presence of a gap for integer spin 
chains, £,(T) remains finite as T ^ 0. We estimate that 
C(l,0)/a = 6.0(1) and ^(2,0)/a = 50(1). These results 
satisfy the relation v = AS, with the values for v and A 
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Fig. 3. The temperature dependence of the staggered sus- 
ceptibility per spin. The solid lines are the results of fits 
to Eq. (8). 



obtained above; As=i/J = 0.40(1), vs=i/Ja = 2.5(1), 
As=.ilJ = 0.090(5), and vs=2/Ja = 4.50(33), and also 
agree with those of previous numerical studies within the 
combined errors [30,34]. For the 5=1/2 chain, we also 
plot in Fig. 2 the WZNW model prediction. 



C(i/2,r) 



ln(To/T) 



ln(ln(To/T) + 1/2) 
21n(ro/T)2 



(7) 

with the parameters obtained by Nomura and Yamada in 
a thermal Bethe ansatz study [35] . 

The staggered susceptibility per spin is shown in Fig. 
3. The different behaviors of integer and half odd integer 
spin chains are clearly manifest in this plot. In Fig. 4, the 
static structure factor Cq at q = tt is plotted as a function 
of temperature. Note, that for integer spin chains (5, T) 
peaks at T « A/ 2: closely similar behavior is observed 
for spin- 1/2 ladders of even width [25]. The extrapolated 
zero-temperature values for the 5=1 chain, Xs(l,0) = 
18.6(1) and C^(1,0) = 3.83(2) agree well with the cor- 
responding results of the exact diagonalization study by 
Sakai and Takahashi [36]. We also obtained Xs(2,G) = 
1160(10), and C^(2,0) = 52.0(3) for the 5 = 2 chain. 

Recently, Starykh, Sandvik, and Singh [37] studied the 
static structure factor and the staggered susceptibility of 
the 5 = 1/2 chain and obtained low-temperature analytic 
forms for these quantities: 

Xs(5,r) = D^iS)T-' [ln(T^(5)/r)]^/^ , (8) 

C^{S, T) = D,{S) [ln{T,{S)/T)f^ . (9) 
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Fig. 4. Static structure factor peak intensity C.,r{S,T) is 

plotted to show the linear dependence of C^^^ on logT. 
The solid lines are fits to Eq. (9) . The arrows indicate the 
peak positions as discussed in the text. 



We are able to fit our 5=1/2 data to these forms and 

thereby to extract D^{\I2) = 0.30(1), T^(l/2) = 9.8(1.2), 
^(1/2) = 0.091(1), and T^{l/2) = 21(1). The solid lines 
in Figs. 3 and 4 are the results of fits to Eqs. (8) and (9) 
for T/J < 0.25. The parameters thus obtained agree with 
those of Ref. [37] to within 5%, except for T^{l/2). This is 
due to the difference in fitting range, since only very low 
temperature data show asymptotic behavior for Xs{S, T). 
Our fitting range is 0.009 < T/J < 0.25, while 0.035 < 
T/J < 0.25 was used in Ref. [37]. Our 5 = 3/2 data 
can also be fitted to Eqs. (8) and (9) for T/J < 0.1. The 
parameters so-obtained are D^{3/2) — 7.8(3), 7^(3/2) = 
3.2(1.1), Ds(3/2) = 0.63(2), and Ts{3/2) = 9.5(1.2). 

As shown in Fig. 2, the correlation length in the spin- 
5/2 system TMMC is described well by the exact clas- 
sical result. From Eq. (7) one sees that quantum effects 
only modify the classical behavior ^ ~ 1/T by a logarith- 
mic correction factor. On the other hand, for 5 = 1/2 
the structure factor in Eq. (9) changes from the classical 
form C^(5,T) - 1/T to the quantum form C^(5,T) ~ 
[ln(Ts(5)/T)]^/^; this is a qualitative and not just a quan- 
titative change due to quantum fluctuations in the di- 
vergence of C7r(5, T) with decreasing temperature. As- 
suming this also holds for 5 = 3/2 and 5 = 5/2, we 
plot C^(5, T)~^ versus T in Fig. 5. It is evident that 
for 5 = 5/2 there is a crossover from classical to quan- 
tum behavior around T/J ^ 1. Interestingly, the exper- 
imental data for TMMC do not seem to exhibit such a 
crossover. However, in TMMC there should also be a spin- 
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Fig. 5. Inverse static structure factor peak intensity 
C7r(<S', T)^^, multiplied by S for graphical purposes, ver- 
sus T/J. Note the deviation of the quantum Monte Carlo 
results for 5 = 5/2 from the data for TMMC around 
T/J~ 1. 



space crossover from Heisenberg to XY behavior at the 
XY gap temperature, which is A^y / J ~ 0.7. In the XY 
regime, €^^{3, T) is enhanced. We speculate therefore that 
in TMMC the quantum and XY effects fortuitously can- 
cel thus leading to the apparent classical 1/T behavior 
observed down to very low temperatures. Future quan- 
tum Monte Carlo calculations for S = 5/2 including the 
XY anisotropy should serve to test this conjecture. 



5 Discussion 

In a recent theoretical study of gapped spin chains at 
non-zero temperature, Damle and Sachdev[27] developed 
a semiclassical picture based on thermally excited parti- 
cles for T <§; /i. In this temperature regime, they were 
able to obtain expressions for several dynamic properties 
of the 0(3) nonlinear a model. They also developed a 
rather different semiclassical approach for A < T < J 
which is based upon classical waves described by the con- 
tinuum 0(3) a model. Using this approach, Damle and 
Sachdev derive one-loop expressions for the uniform sus- 
ceptibility and correlation length in the temperature range 
A<T<J: 
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Fig. 6. (a) Correlation length and (b) uniform suscepti- 
bility as a function of temperature for integer spin chains. 
The solid lines are the results Eqs. (10) and (11) without 
any adjustable parameter. The dashed line in (a) is the 
relation v = A^, and the dashed line in (b) is Eq. (4). 



vs 
2-kT 



( Z2'kT 



(11) 



where 7 = 0.5772... is Euler's constant. 

In Fig. 6, we show T^{S.T)/vs and vsXu{S,T) as 
functions of T /As without any adjustable parameters, 
since As and vs have been determined independently. 
For the S* = 1 chain we use values As=i = 0.41050(2), 
C(l,0)/a = 6.03(1), and vs=i/Ja = 2.49(1) taken from 
the literaturc[29,30], since these zero temperature results 
have smaller error bars than our own values. The solid 
lines are Eqs. (10) and (11), which agree well with the 
S — 2 data for As=2 < T < J. An even better agreement 
might be obtained if the theory were extended beyond the 
one loop approximation. It is evident, nonetheless, that a 
window of temperature in which Eqs. (10) and (11) apply 
indeed exists for the S = 2 chain, as speculated by Damle 
and Sachdev [27]. 

In Fig. 6(b), we also show the scaling of the uniform 
susceptibility at low temperatures (T <C A) along with the 
theoretical expression Eq. (4) as a dashed line without any 
adjustable parameters. The dashed line in Fig. 6(a) shows 
that the relation v = A^ holds up to T Zi/2 in integer 
spin chains. 

Spin- 1/2 ladders are arrays of n coupled Heisenberg 
chains with S = 1/2. In analogy to integer spin chains, 
ladders with an even number n of chains exhibit exponen- 
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Fig. 7. Uniform susceptibility per spin and staggered 
correlation length of the quantum spin chains are com- 
pared with those of quantum spin ladders, consisting 
of n isotropically coupled antiferromagnetic S = 1/2 
chains[25]: (a) correlation length of spin chains; (b) corre- 
lation length of spin ladders; (c) uniform susceptibility of 
spin chains; and (d) uniform susceptibility of spin ladders. 
Solid lines in (a) and (b) are the WZNW model prediction, 
Eq. 7 

tially decaying correlations in their ground state due to 
the presence of a spin gap, while those with odd n show 
behavior similar to that of half-odd-integcr spin chains. 
Therefore, it is illuminating to plot the results such that 
the behaviors of quantum spin chains can be qualitatively 
compared with those of 5 = 1/2 quantum spin ladders. 
In Fig. 7, we display the staggered correlation length and 
the uniform susceptibility per spin for both spin chains 
and isotropically coupled spin ladders[25]. The same sym- 
bol is used for the spin-S" chain and the spin-1/2 ladder of 
width n = 2S. One observes the markedly different behav- 
iors of gapped systems shown in open symbols compared 
with those of gapless systems shown as solid symbols. 

In summary, antiferromagnetic Heisenberg spin chains 
with spin values ranging from 5=1/2 to 5 = 5/2 have 
been studied with the quantum Monte Carlo method. The 
temperature dependences of the uniform susceptibility, 
the staggered susceptibility, the static structure factor peak 
intensity, and the correlation length are obtained. We find 
that at high temperatures these quantities agree very well 
with the exact results for the classical spin chain, and that 
quantum effects become progressively more important as 



the temperature is decreased. In addition, our data for Xu 
and ^ for the 5 = 2 chain in the intermediate temperature 
range As=2 <T < J are reasonably well predicted by the 
theory of Damle and Sachdev. 
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